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Abstract 

This paper demonstrates how well the k-cn turbulence 
model aescribes the nonlinear growth of flow instabilities 
from laminar flow into the turbulent flow regime. Viscous 
modifications are proposed for the k-cj model that yield 
close agreement witn measurements and with Direct 
Numerical Simulation results for channel and pipe flow. 
These modifications permit prediction of subtle sublayer 
details such as maximum dissipation at the surface, k ^ y 2 
as y 0, and the sharp peak value of k near the surface. 
With two transition specific closure coefficients, the model 
equations accurately predict transition for an incompressi- 
ble flat-plate boundary layer. The analysis also shows why 
the k-c model is so difficult to use for predicting transition. 


I. Introduction 


mass transfer. The purpose of this paper is to build upon 
the work of Wilcox, et al, taking advantage of recent Direct 
Numerical Simulation (DNS) results in developing appro- 
priate viscous modifications for the Wilcox 7 k-u model. 
The scope of this paper is confined to incompressible flows. 


For incompressible turbulent fluid flow, the complete set 
of equations that constitute the Wilcox k-co two-equation 
model are as follows. 
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There has been renewed interest in development of 
methods for predicting boundary-layer transition. Current 
interest in vehicles such as the National Aerospace Plane 
(NASP), for example, provides the impetus for developing 
accurate transition prediction tools. Furthermore, because 
hypersonic boundary layers rarely achieve momentum- 
thickness Reynolds numbers large enough to sustain fully- 
developed turbulence, even the post-transition region 
generally exhibits nontrivial viscous effects. Consequently, 
accurate low-ReynoIds-number turbulence models are also 
needed. 

The standard approach is to view development of a 
transition model and a low-Reynolds-number turbulence 
model as two separate issues. The strongest argument in 
favor of this approach is simply that all spectral effects are 
lost in the time-averaging process used by turbulence 
models. Tollmien-Schlichting waves, for example, cannot 
be distinguished by a turbulence-model. Since a given 
boundary layer is unstable to perturbations that fall in a 
specific range of frequencies, conventional turbulence 
models, which distinguish only magnitude and an average 
frequency, can never be certain if a given perturbation will 
actually cause transition. However, if we implement two 
separate models, one for the transition region and another 
for the developing turbulent region, achieving a smooth 
joining of the two models’ predictions presents an addition- 
al complication. 

This complication can be avoided if we view both issues 
as low-Reynolds-number phenomena that can be addressed 
in the context of a single model. The strongest argument 
for this approach is that we can use the same model to 
describe a smooth transition from laminar to fully turbulent 
flow, including the transitional region. This approach is 

E lausible provided we restrict our applications to broad- 
and transition-triggering disturbances. 

The research of Wilcox, et al 1 ^ provides a great deal of 
support for the latter approach. Using a k-oS 2 turbulence 
model and transition-specific, low-Reynolds-number modi- 
fications, Wilcox simulated boundary-layer transition for a 
wide range of Mach numbers including pressure gradient, 
surface roughness, surface heating and cooling, and surface 
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In Equations (1-8), t is time, x, is position vector, u, is 
velocity, p is density, p is pressure, u is molecular viscosity, 
and ^ is the sum o? molecular and Reynolds stress tensors. 
Also, £jj is the Kronecker delta, k is the turbulence kinetic 
energy, is specific dissipation rate, is Reynolds stress 
tensor, and ^ is eddy viscosity. The six parameters a, a*, 
0, 8*, a and o* are closure coefficients whose values are 
given below. 
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where Re T is turbulence Reynolds number defined by 
Re T = p k / (gj ^) (15) 

Section II explains in detail how the turbulence model 
simulates transition, and justifies the form of the viscous 
modifications. Section III explains why the most popular 
low-Reynolds-number k-e models are unsuitable tor pre- 
dicting transition. Section IV demonstrates how well the 
model performs for low-Reynolds-number channel flow 
and for pipe flow. Section V includes transition predictions 
for an incompressible flat-plate boundary layer. Section VI 
presents a summary of and conclusions drawn from the 
study. 


II. Simulatin; 


ig Transition yyith Turbulence 
Model Equations 


Turbulence model equations can be used to predict 
transition from laminar to turbulent flow, although most 
models predict transition to turbulence at Reynolds num- 
bers that are at least an order of magnitude too low. To 
understand why and how the k-co model predicts transition, 
consider the flat-plate boundary layer. For the k-oi model, 
the incompressible, two-dimensional boundary-layer form 
of the equations for k and co is as follows. 
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There are two important observations worthy of mention 
at this point. First, if the turbulence energy is zero. Equa- 
tion (21) has a well-behaved solution. That is, when k = 0, 
the eddy viscosity vanishes and the co equation uncouples 
from the k equation. Consequently, the k-oo model has a 
nontrivial laminar-flow solution for co. Second , the signs of 
P|c and determine whether k and co are amplified or re- 
duced in magnitude. However, it is not obvious by inspec- 
tion of Equations (22) and (23) how the signs of these terms 
vary with Reynolds number as we move from the plate 
leading edge to points downstream. We can make the 
variation obvious by rewriting Equations (22) and (23) in 
terms of the Blasius transformation. 

Before we introduce the Blasius transformation, we must 
determine the appropriate scaling for c o. To do this, we 
note that close to the surface of a flat plate boundary' layer, 
the specific dissipation rate behaves according to 7 
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In terms of the Blasius similarity variable, t), defined by 

n = y/(vx/U.)^ (25) 

where U w is freestream velocity, the asymptotic behavior of 
co approaching the surface is 
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Consequently, we conclude that the appropriate scaling for 
co in the Blasius boundary layer is given by 


co = W( x, ) 
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where W(x,iri is a dimensionless function to be determined 
as part of the solution. Hence, if we write the velocity in 
terms of dimensionless velocity, F(x,tj), i.e., 

u - U. F(x,n) (28) 

the net production per unit dissipation terms become 


where u and v are velocity components in the streamwise 
(x) and normal (y) directions, respectively, v is kinematic 
molecular viscosity, and v T is kinematic eddy viscosity. We 
can most clearly illustrate how the model equations predict 
transition by rearranging terms in Equations (17) and (18) 
as follows. 
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The net production per unit dissipation for the two 
equations, Pfc and P^, are defined by: 
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Thus, both P^ and P Q increase linearly with Reynolds 
number. Re,. From the exact laminar solution for F(tj) and 
W(n) [the x dependence vanishes for the Blasius boundary 
layer], the maximum value of the ratio of 3F/ 3n to W is 
given by 
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Hence, as long as the eddy viscosity remains small com- 
pared to the molecular viscosity, we can specify the precise 
points where P k and change sign. In general, using 
Equation (31), we conclude that the sign changes occur at 
the following Reynolds numbers. 


objective is to match the minimum critical Reynolds 
number. Reference to Equation (32) shows that we must 
require 

B*/a* -► 1 as Re T -► 0 (35) 
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With no viscous modifications, the closure coefficients a, 
a*, & and 0* are 5/9, 1, 3/40 and 9/100, respectively. 
These values correspond to the limiting form of Equations 
(9-14) as Re x -* ». Using these fully turbulent values, we 
find (Re x ) k * 8,100 and (Re = 12,150. Thus, starting 
from laminar flow at the leading edge of a flat plate, the 
following sequence of events occurs. 


(1) The computation starts in a laminar region with 
k = 0 in the boundary layer and a small freestream 
value of k. 

(2) Initially, because P k < 0 and P w .< 0, dissipation of 
both k and co exceeds production. Turbulence 
energy is entrained from the freestream and spreads 
through the boundary layer by molecular diffusion. 
Neither k nor co is amplified and the boundary layer 
remains laminar. 

(3) At the critical Reynolds number , Rev,. = 8,100, 
production catches dissipation in the k equation. 
Downstream of k production exceeds k dissipa- 
tion and turbulence energy is amplified. At some 
point in this process, the eddy viscosity grows rapidly 
and this corresponds to the transition point. 

(4) k continues to be amplified and, beyond 
Re x = 12,150 production catches dissipation in the w 
equation. u> is now amplified and continues growing 
until a balance between production and dissipation 
is achieved in the k equation. When this balance is 
achieved, transition from laminar to turbulent flow 
is complete. 


Consistent with experimental measurements, the entire 
process is very sensitive to the freestream value of k. There 
is also a sensitivity to the freestream value of go, although 
the sensitivity is more difficult to quantify. 

Three key points are immediately obvious. First , k 
begins growing at a Reynolds number of 8,100. By contrast, 
linear-stability theory tells us that Tollmien-Schlichting 
waves begin forming in the Blasius boundary layer at a 
Reynolds number of 90,000. This is known as the mini- 
mum critical Reynolds number . Correspondingly, we 
find that the model predicts transition at much too low a 
Reynolds number. Second , inspection of Equations (32) 
and (33) shows that the width of the transition region is 
controlled by the ratio of 6 to aa*. Third . transition will 
never occur if reaches zero earlier than P k . Thus, 
occurrence of transition requires 

aa * < a m B/B* as Re T -*■ 0 (34) 

This fact must be preserved in any viscous modification 
to the model. The viscous modifications in Equations (9- 
14), i.e„ the dependence of a, a* and 0* upon Re x , are de- 
signed to accomplish two objectives. The most important 


A secondary objective is to achieve asymptotic consist- 
ency with the exact behavior of k and dissipation, 
e = 0*kco, approaching a solid boundary. That is, we would 
like to have 


k/y 2 constant, e/k -* 2v/y 2 as y -* 0 (36) 

Close to a solid boundary, Wilcox 7 shows that the dissipa- 
tion and molecular diffusion terms balance in both the k 
and 0 ) equations. The very-near-wall solution for is given 
by Equation (24). The solution for k is of the form 

k/y n -► constant as y 0 (37) 

where n is given by 

n = Vz[\ + (1 + 240* /0)^) (38) 

Noting that dissipation is related to k and to by 
e = 0*ko> (39) 


we can achieve the desired asymptotic behavior of Equa- 
tion (36) by requiring 

0*/B - 1/3 as Re x - 0 (40) 

Requiring this limiting behavior as Re x 0 is sufficient to 
achieve the desired asymptotic behavior as y 0 since the 
eddy viscosity, and hence, Re^- vanish at a solid boundary. 


If we choose to have 6 constant for all values of Re T , 
Equations (34), (35) and (40) are sufficient to determine 
the limiting values of a* and 0* and an upper bound for 
aa* as turbulence Reynolds number becomes vanishingly 
small. Specifically, we find 


aa * < 0 

a* -* 0/3 
0 * - 0/3 


as Re x — 0 (41) 


Wilcox, et al, M make the equivalent of aa * and a* in 
their k-GJ 2 models approach the same limiting value and 
obtain excellent agreement with measured transition width 
for incompressible boundary layers. Numerical experimen- 
tation with the k-c*) model indicates the optimum choice for 
incompressible boundary layers is aa * 0.74/8, or 


aa * -*> 1/18 as Re x -*■ 0 


(42) 


Equations (9-14) postulate functional dependencies 
upon Re x that guarantee the limiting values in Equations 
(41) and (42), as well as the original fully turbulent values 
for Re x -» x. 


The three coefficients R e , R k and R^ control the rate at 
which the closure coefficients approach their fully turbulent 
values. As in previous analyses based on the k-o) model, 7 * 89 
we can determine their values by using perturbation 
methods to analyze the viscous sublayer. Using the well 
established procedure, we can solve for the constant in the 
law of the wall, B, by solving the sublayer equations and 
evaluating the following limit. 

r 1 

B = lim [u + - - Jin y + ] 
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where u + = u/i^ and y + = i^y/v are standard sublayer 
scaled coordinates. Also, k = 0,41 is Karman’s constant. 
For a given value of R a and R k , there is a unique value of 
R w that yields a constant in the law of the wall of 5.0. For 
small values of R s the peak value of k near the surface is 
close to the value achieved without viscous corrections, viz, 
u^/ffl*)*. As R s increases, the maximum value of k near 
the surface increases. Comparison of computed sublayer 
structure with Direct Numerical Simulation (DNS) results 
of Mansour, et al 10 indicates the optimum choice for these 
three coefficients is as indica*ed in Equation (14). Section 
IV presents a complete comparison of computed channel 
flow properties with the Mansour, et al DNS results. 

The only flaw in the model’s asymptotic consistency 
occurs in the Reynolds shear stress, t^. While the exact 
asymptotic behavior is 'v y 3 , the model as formulated 
predicts r ^ y 4 . This discrepancy could easily be removed 
with another viscous modification. However, results ob- 
tained to date indicate this is of no significant consequence. 
It has no obvious bearing on either the model’s ability to 
predict transition or properties of interest in turbulent 
boundary layers. The additional complexity and uncertain- 
ty involved in achieving this subtle feature of the very-near- 
waU behavior of does not appear to be justified. 


Equation (48) is nearly identical to the limiting form of 
Equation (18) for v T /v -► 0. The only significant difference 
is the last term on the right-hand side of Equation (48). 
Except dose to the surface where k must be exactly zero, 
this term is unlikely to have a significant effect on the solu- 
tion for small nonzero values of k. However, in a numerical 
solution, products of dependent-variable gradients are 
generally destabilizing, ana the problem can only be aggra- 
vated by having a coeffident inversely proportional to k. 
This is not an insurmountable problem. However, estab- 
lishing starting conditions is dearly more difficult with the 
k-f model than with the k-<u model. 

Given the diverse nature of viscous modifications that 
have been proposed for the k-e model 12 it is impossible to 
make any universal statements about why a specific model 
fails to predict realistic transition Reynolds numbers. 
Perhaps the strongest statement that can be made is, no one 
has approached the problem from the transition point 
of view. Most researchers have sought only to achieve 
asymptotic consistency and attempted transition predictions 
only as an afterthought. We can gain some insight by 
examining the net production per unit dissipation terms for 
the k and € equations that are analogous to Equations (29) 
and (30), viz, 


III. Difficulties Attending Use of th e k-f 

Model 

Given the information developed in Section II, it is a 
simple matter to explain why little progress has been made 
in predicting transition with the k-e model. 11 The primary 
difficulties can be easily demonstrated by focusing upon 
incompressible boundary layers. If we use the standard 
form of the k-€ model. Equations (17-19) are replaced by 
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Values of the closure coefficients/functions C fl , C f j, C K . 
f a k and o ( differ amongst t he various versions of this 
model. In the absence of viscous modifications, the stand- 
ard values used are C fl = 1.44, C e2 = 1-92, C M = .09, f^ = 1, 
a k = 1 and o t = .769. 

One critical difference from the k-co model is obvious by 
inspection of Equations (44-46). Specifically, if the turbu- 
lence energy is zero, c must also be zero. We cannot simply 
drop the eddv viscosity in the f equation because^ of the 
presence of k in the denominator or the e equation’s dissi- 
pation term. The model does possess a laminar-flow solu- 
tion for the ratio of e to k. That is, if we make the formal 
change of variables 

e = C„kco 

and assume v T < < v, the following laminar-flow equation 
for results. 
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On the one hand , without viscous damping, if we 
assume Equation (31) is valid, we find (Re x )|^ = 8,100 and 
(Re x ) f = 10,800. Consequently, like the high-Reynolds- 
number version of the k-oo model, transition will occur at 
too low a Reynolds number. 

On the other hand , because f^, C (2 and sometimes C ' j 
are often permitted to be functions of distance from the 
surface and/or functions of Re x , we cannot simply use 
Equation (31). Furthermore, some modelers add t erms to 
the k and e equations in addition to damping the jtosure 
coefficients. Each set of values for the closure coefficients 
and additional terms must be used in solving Equation (48) 
to determine the laminar-flow solution for e/k. While it is 
cleariv impossible to make a quantitative evaluation of all 
variants of the k-e model, we can nevertheless make some 
general observations. 

From the analysis of the k-oi model, it is obyjqusjhat 
having f < 1 will tend to delay transition. Virtually all 
modelers implement an f„ that will accomplish this end. 
However, the modifications of Jones and Launder, 11 
Chien, 13 and Lam and Bremhorst, 14 for example, damp C (2 
to the extent that (Re x ) e is smaller than (Re x )k. This is the 
opposite of what is needed and will have an undesirable 
effect on both the onset of and the extent of the transition 
region. 

This discussion is not intended as an exhaustive survey of 
the numerous low-Reynolds-number versions of the k-6 
model. Rather, it is intended to illustrate how difficult it is 
to apply the model to the transition problem. Given 
enough additional closure coefficients and damping func- 
tions, the k-6 model can probably be modified to permit 
satisfactory transition predictions. However, even if this is 
done, establishing starting conditions will ultimately require 
a solution to Equation (48). That is, to initialize the compu- 
tation, we must effectively transform to the k-to model. 
Since this is the natural starting point, it seems illogical to 
perform subsequent computations in terms of k and e. 
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IV. Turbulent Flow Applications 

To achieve a complete description of the transition from 
laminar to turbulent flow, we must be able to accurately 
describe the flow in the turbulent regime. This is, after all, 
the primary advantage of using turbulence model equations 
to describe transition. In this section, we examine channel 
and pipe flow to demonstrate how well the low-Reynoids- 
number form of the k-u) model predicts properties or turbu- 
lent flows. 

Figures 1 and 2 compare computed channel-flow skin 
friction, c f with the Halleen and Johnston 15 correlation for 
Reynolds number based on channel height, H, and average 
vefocity ranging from 1G 3 to 10 5 . The correlation is 

Cf = 0.0706 Re H ’^ 4 (51) 

As shown, computed c f differs from the correlation by less 
than 3 % except at the lowest Reynolds number shown 
where the correlation probably is inaccurate. Figure 3 
compares computed pipe flow c f with Prandtl’s universal 
law of friction, viz, 

c { -Vi = 4 Iog 10 (2 Re D c/ 2 ) - 1.6 (52) 

Reynolds number based on pipe diameter, D, and average 
velocity varies from 10 3 to 10 6 . As with channel flow, 
computed c f falls within 59c of the correlation except at the 
lowest Reynolds number indicated where the correlation is 
likely to be in enor. 

For more detailed comparisons, we consider two low- 
Reynolds-number channel-flow cases corresponding to the 
DNS results of Mansour, et al 10 and one nigh-Reynolds- 
number pipe flow case corresponding to measurements of 
Laufer 16 . 


For purposes of identification, the three cases are re- 
ferred to in terms of the parameter 

Re r = u, R / v (53) 

where u* is friction velocity and R is either channel half 
height or pipe radius. Figures 1 and 2 compare various 
computed profiles with the Mansour, et al DNS results for 
Re f = 180 and 395, respectively. 

Six different comparisons are shown in each figure, 
including mean velocity, skin friction. Reynolds shear 
stress, turbulence kinetic energy, turbulence energy produc- 
tion and dissipation rate. For both cases, velocity, Reynolds 
shear stress, and turbulence kinetic energy profiles differ by 
less than 7%. Most notably, for both Reynolds numbers, 
the model predicts the peak value of k near the channel 
wall to within 4% of trie DNS value. Additionally, ap- 
proaching the surface, the turbulence-energy production, 
r^SU/Oy, and dissipation, e, are within 10% of the DNS 
results except very dose to the surface. 

Figure 3 compares k-cu model pipe flow results with 
Lauferis 16 measurements at a Reynolds number based on 
pipe diameter and average velodty of 40,000. As shown, 
computed and measured velodty and Reynolds shear stress 
profiles differ by less than 89c. As with channel flow, 
computed and measured turbulence kinetic energy differ by 
about 59c including close to the surface where the sharp 

r iak occurs. Note diat, at this high a Reynolds number, the 
profile has a sharp spike near y = 0 and this feature is 
captured in the computations. Except very close to the 
surface, computed turbulence energy production and dissi- 
pation differ from measured values by less than 10%. This 
may actually be a desirable result. That is, some controver- 
sy exists about the accuracy of Laufers dissipation meas- 
urements dose to the surface. 



Figure 1. Incompressible channel flow with Re f = 180, Figure 2. Incompressible channel flow with Re T = 395, 
Re H s 5,590; o Mansour, et al; ■ Johnston Formula. Re H = 13,750; o Mansour, et al; ■ Johnston Formula. 
















Figure 3. Incompressible pipe flow with Re r = 1058, 
Re H = 40 t 000; o Laufer; ■ Prandtl Formula. 


V. Transition Applications 

Figure 4 compares computed and measured transition 
Reynolds number, Re$ t , tor an incompressible flat-plate 
boundary layer. We define the transition Reynolds number 
as the point where the skin friction achieves its minimum 
value. Results are displayed as a function of freestream 
turbulence intensity, T, defined by 

Y " ,0O( fc7 ) * (54) 

where subscript e denotes the value at the boundary-layer 
edge. As shown, consistent with the data compiled by 
Dryden, 17 Re$ ( increases as the freestream intensity de- 
creases, Because cu can be thought of as an averaged 
frequency of the freestream turbulence, it is reasonable to 
expect the predictions to be sensitive to the freestream 
value of co. To assess the effect, the freestream value of the 
turbulence length scale defined by Jt = k*vcu has been 
varied from .0015 to .1005 where 5 is boundary layer thick- 
ness. As shown, computed Re e t values bracket virtually all 
of the data. These predictions are markedly superior to the 
preliminary efforts of Wilcox 18 in developing iow-Reynolds- 
nuinber modifications for the k-co model. 

Figure 5 compares computed width of the transition 
region with measurements of Dhawan and Narasimha. 19 
We define transition width, Ax t , as the distance between 
minimum and maximum skin-friction points. The comput- 
ed width fails within experimental data scatter for 
Re Xt < 10 6 , and lies a bit above the data for larger values. 
Ax t is unaffected by the freestream value of w. 



Figure 4. Comparison of computed and measured variation 
of transition Reynolds number with freestream turbulence 
intensity; incompressible flat-plate boundary layer. 





Figure 5. Comparison of computed and measured width of 
the transition region for an incompressible flat-plate 
boundary layer. 


VI. Summ ary and Conclusions 

The primary objective of this paper has been to illustrate 
how two-equation turbulence models, most notably the k-oj 
model, predict transition. While the long-term goal of this 
research is to develop a transition model for all Mach 
numbers, this paper has focused on the case about which we 
know the most, viz, incompressible flow. The low- Reyn - 
olds-number modifications proposed not only facilitate 
accurate transition predictions, but also yield reasonably 
close agreement with DNS results for low-Reynolds- 
number channel flow. 

The degeneracy of the € equation in the k-e model is a 
major stumbling block that impedes successful application 
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to the transition problem. By transforming the model to an 
equivalent k-oo model, it would be possible to remove some 
of the difficulties. After making such a transformation 
however, there is little reason to transform back. 

The applications presented in Sections IV and V indicate 
we have formed a solid foundation for future low-Reyn- 
olds-number and transition research. 


VII- Acknowledgements 

Research was supported by the U. S. Army Research 
Office and the NASA Langley Research (Tenter under 
Contract DAAL03-89-C-0032 with Dr. Thomas Doligalski, 
Dr. Julius Harris and Dr. Thomas Zang as Contract Moni- 
tors. 


References 

Wilcox, D. C„ Turbulence-Model Transition Pre- 
dictions, ” AIAA Journal, Vol. 13, No. 2, pp. 241-243 (1975). 

2 Wilcox, D. C, 'Turbulence Model Transition Pre- 
dictions: Effects of Surface Roughness and Pressure Gradi- 
ent,*' AIAA Paper 75-857, Hartford, CT (1975). 

Wilcox, D. C. and Chambers, T. L„ "Application of 
the Turbulence-Model Transition-Prediction Method to 
Flight Test Vehicles," Turbulence in internal Fl o w s. S.N.B. 
Murthy, pp. 233-247 (1976). 

4 Wilcox, D. C., "A Model for Transitional Flows," 
AIAA Paper 77-126, Los Angeles, CA ( 1977). 

Wilcox, D. C., "Alternative to the e 9 Procedure for 
Predicting Boundary-Layer Transition," AIAA Journal, Vol. 
19, No. 1, pp. 56-64 (1981). 

Wilcox, D. C. and Marvin, J. G., "Combined Effects 
of Freestream Turbulence and Mass Addition on Blunt- 
Body Heating and Transition," Proceedings of the 2 nd 
ASME/JSME Joint Thermal Engineering Conference, 
Honolulu, HI (1987). 

Wilcox, D. C., "Reassessment of the Scale Deter- 
mining Equation for Advanced Turbulence Models,” AIAA 
Journal, Vol. 26, No. 11, pp. 1299-1310(1988). 

Wilcox, D. C. and Traci, R. M., "A Complete Model 
of Turbulence," AIAA Paper 76-351, San Diego, CA (1976). 

Wilcox, D. C. and Rubesin, M. W., "Progress in 
Turbulence Modeling for Complex Flow Fields Including 
Effects of Compressibility," NASA TP 1517 (1980). 

10 Mansour, N. N., Kim, J. and Moin, P., "Reynolds 
Stress and Dissipation Rate Budgets in Turbulent Channel 
Flow," JFM, Vol. 194, pp. 15-44 (1988). 

1 Hones, W. P. and Launder, B. E., ’The Prediction of 
Laminarization with a Two-Equation Model of 
Turbulence,” International Journal of Heat and Mass 
Transfer, Vol. 15, pp. 301-314 (1972). 

l2 Patel, V. C., Rodi, W. and Scheuerer, G., Turbu- 
lence Models for Near-Wall and Low Reynolds Number 
Flows: A Review," AIAA Journal, Vol. 23, pp. 1308-1319 
(1985). 


14 Lam, C. K. G. and Bremhorst, K. A., "Modified 
Form of the k-e Model for Predicting Wall Turbulence 
Journal of Fluids Engineering, Vol. 103, pp. 456-460 (1981). 

l5 Halleen. R. M. and Johnston, J. P., The Influence 
of Rotation on Flow in a Long Rectangular Channel - An 
Experimental Study/’ Report MD-18. Mechanical Engi- 
neering Dept., Stanford University (1967). 

16 Laufer, J., The Structure of Turbulence in Fully 
Developed Pipe Flow," NACA 1174 (1952). 

17 Dryden, H. L., Aerodynam ics and Jet Propulsion, 
Vol. V, University Press, Princeton, NJ (1959). 

iWilcox, D. C., "Progress in Hypersonic Turbulence 
Modeling," AIAA Paper 91-1785, Honolulu, HI ( 1991). 

19 Dhawan. S. and Narasimha, R., "Some Properties 
of Boundary Layer Flow During the Transition from 
Laminar to turbulent Motion," JFM, Vol. 3, pp. 418-436 
(1958). 


13 Chien, K-Y., "Predictions of Channel and Bound- 
ary-Laver Flows with a Low-Reynolds-N umber Turbulence 
Model," AIAA Journal, Vol, 20, pp. 33-38 ( 1982). 





